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Abstract 

We study the Jahn-Teller effect of positive fullerene ions 2 C^ and l C$. The 
aim is to discover if this case, in analogy with the negative ion, possesses a 
Berry phase or not, and what are the consequences on dynamical Jahn-Teller 
quantization. Working in the linear and spherical approximation, we find no 
Berry phase in 1 Cg ( |, and presence/absence of Berry phase for coupling of one 
L = 2 hole to an L = 4/L = 2 vibration. We study in particular the special 
equal-coupling case (#2 = 9i), which is reduced to the motion of a particle 
on a 5-dimensional sphere. In the icosahedral molecule, the final outcome 
assesses the presence/absence of a Berry phase of ir for the h u hole coupled to 
Gg/Hh vibrations. Some qualitative consequences on ground-state symmetry, 
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low- lying excitations, and electron emission from Ceo are spelled out. 
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I. INTRODUCTION 



Non-adiabatic systems, where the quantum mechanics of electrons and that of ions are 
deeply entangled, are among the most fascinating in physics. In condensed matter, entire 
classes of phenomena such as superconductivity or charge- density waves are of this nature. 
Dynamical Jahn- Teller (DJT) centers (molecules and ions, impurities, nuclei) constitute per- 
haps the simplest "zero-dimensional" example of nonadiabatic behavior, and have attracted 
attention for over half a century M. 

An amusing anomaly has long been pointed out in certain simple DJT centers, where 
an e electronic doublet is Jahn- Teller (JT) coupled to a vibrational ( "pseudo- rotational" ) E 
doublet. In modern language, this anomaly can be seen as a Berry phase fHH. In simple 
words, what happens in these systems is that there exists a class of closed paths (loops) 
in the E coordinate space, regarded as classical, adiabatically transporting the electronic 
ground state from an initial value ipo to a final value e 17 ^ — ~~ V'oj instead of into ipo itself. 
There is a Berry phase 7 = 7r, whose origin lies in the fact that the point of electronic 
degeneracy (the "conical point" [7|,||) is irreducibly surrounded by the trajectory in vibron 
space. 

This anomaly which is present at the adiabatic level, where the ionic coordinated are 
classical, has very interesting consequences when the full joint ion-electron problem is finally 
quantized. In the e <g) E case cited above, there appears a new conserved "pseudo-angular 
momentum" quantum number, whose values are half- integer Jl]|| . Consequently, all vibronic 
levels become twofold degenerate, as appropriate to a rotor coupled to a pseudospin | . 

It is an interesting question therefore, to inquire whether other nontrivial cases of DJT are 
realized in nature, which also possess a Berry phase. High-symmetry polyatomic molecules 
are an obvious place to look at. 

The fullerene molecule, Ceo, has a very high X/j symmetry, and a nondegenerate ground 
state. Because of high symmetry, however, fullerene ions, both negative and positive have a 
high orbital degeneracy, and are therefore good candidates. 
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In our previous work, we have studied the DJT problem of negatively- charged 
fullerene ions, where 1 < n < 5 electrons partly fill the t lu Lowest Unoccupied Molecular 
Orbital (LUMO) of the neutral molecule. This electronic state couples linearly to the five- 
fold degenerate H g vibrations. For odd electron number n, the spectrum of this coupled 
dynamical t lu ® H g system is affected by anomalies which can again be usefully described 
in terms of a Berry phase J5J. 

Among notable consequences, we have predicted the dominance of p-wave attachment 
for low-energy electrons |l0|| , a negative pairing energy for odd n HfJ, characteristic spec- 



troscopic anomalies fllT] and a mechanism for kinematic electron pairing ||12| 

From the mathematical point of view, nevertheless, the DJT of Cg , namely t <S> H, was 
not particularly new: the analytical machinery was basically the same as that previously 
studied by O'Brien for the t <8> (E © T) equal coupling case. 

In the present paper, we move on to analyze the problem of the JT effect in the positively 
charged fullerene ions Cgg", inquiring in particular about the possible relevance of a Berry 
phase in these systems too. As it will turn out, the answer is again affirmative (though with 
some qualifications). Moreover, we shall deal in this case with the less trivial JT problems 
h®H and h®(G@ H), where the search for an adiabatic anomaly will now be done afresh. 

Fullerene cations Cgg" are obtained by removing n < 10 electrons from (i.e. adding 
n < 10 holes to) the fivefold degenerate h u Highest Occupied Molecular Orbital (HOMO) 
|T7Hl9[. Our approach allows us to treat the one-hole case (n = 1) as well as the two-holes 
(n = 2) singlet case, in a way which is essentially analytical, albeit for a very idealized set 
of parameters. 

A hole in the HOMO interacts with the Ceo vibrational modes. We focus on modes which 
couple linearly with the hole density fluctuations (which in turn are of course quadratic in 
the hole wave function). Hence, the icosahedral symmetry of the molecule restricts the 
linearly coupled modes to those of the same symmetry as the irreducible representation of 
the group (Z/J contained in the symmetric product of a hole representation with itself |EU 
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{K ® M, ~ S © © (& © ^) (!) 

The Xft group is a subgroup of the full rotation group 0(3). Therefore, a general represen- 



tation of 0(3) decomposes into a finite number of representations of pl |. In particular, 
the L = 2 £>( 2± ) representations correspond to h g / u in icosahedral symmetry, and L = 4 
X>( 4± ) decomposes into <7 ff / u © h g / u . Equation (|I|) may be rewritten, in 0(3) notation, 

{V^ <g> p( 2 -)} s « © X?( 2+ > © D< 4+ ) (2) 

In order to take advantage of spherical symmetry, and as a simplifying approximation at 
this initial stage, we neglect the icosahedral details of the molecular shape, and we study the 
hole-vibron coupling within the 0(3) scheme, as though Cgo was a perfect sphere |[22|| . We 
study here therefore the problem of T>^ 2 ~" > holes coupled to T>^ 0+ \ £> ( - 2+ - ) and X>( 4 +) vibrations. 

The coupling with a X>( 0+ ) vibration is a trivial polaronic problem, solved exactly as 
a displaced oscillator, and is irrelevant here, since it shifts energies, but does not cause 
splittings. We concentrate on quadrupolar {T>^ 2+s> ) and hexadecapolar (X'^ 4+ - ) ) modes, which 
lead instead to a non-trivial JT effect. In order to simplify the notation, we shall switch 
from the 0(3) to the SO (3) notation, omitting therefore the ± signs for the inversion. 

The Hamiltonian we consider has the following structure: 

H = H + H h - V , Hh- V = H^_ v + H^2 V , (3) 

where Hq describes a free (uncoupled) hole and T>^ 2 \ T>^ vibrations, and Hh- V the 
hole-vibron interaction, namely 

L ( 1\ 2 

H = huJ L ( b L,m b L,m + r ) + (e ~ fl) ^ £ 4n,a C m,a , 



L=2,4 m=-L x ~ cr=t,| m=-2 

T 



H { h L l = 9L^ £ (-ir 2 (^^i|2,-m 2 ;2,m 1+ m 2 )x (4) 



<r,mi ,m2 



Here, 7tu;£ are the energies of the (harmonic) vibrations &^ m , g"£ are dimensionless coupling 
parameters, < L,m | Li, mi; L2, ^2 > are Clebsch-Gordan coefficients, and there is no loss 



of generality in choosing the energy zero so that e = //. For simplicity, we shall assume here 
a single L = 2 and a single L = 4 mode corresponding to £>( 2 )(= and T>^\= G g + H g ) 
vibrons respectively. In real Ceo, the spectrum is more complex (see Tab. |), involving 8 
H g and 6 G g modes. Of these, only two H g modes, namely H g (l) and H g {4) are actually 
derived from a T>^ spherical vibration, while only two doublets, namely G g (l) @H g {2) and 
G g (4:) © H g (6) are of L = 4 origin The other modes derived from different angular 

momenta (L = 3,5,6,7) do not belong to a complete L-multiplet, and cannot be simply 
treated in spherical symmetry. Our work is therefore incomplete and of purely qualitative 
value. Nevertheless, we expect some of the features of the present simplified analysis to be 
of more general relevance. For example, close pairs of modes, such as G g (3) © H g (3), or 
G g (6) © H g (8), have no reason to behave qualitatively differently from our T>^ mode, in 
spite of their different L = 6 and L = 7 origins. We will also find a "conservation of Berry 
phases": if a JT-coupled problem has a Berry phase, the extension of the vibron space 
with the addition of a mode that does not carry a Berry phase, leaves intact the original 
Berry phase. In our case, we will show that the Berry phase is coming exclusively from the 
h u © Gg coupling, the addition of the H g modes changing the dynamical properties, but not 
the topological ones. 

According to the different relative values of the four parameters (u; 2 , ^2, ^4, (74) in H, 
different regimes can be attained. In particular, when the coupling is extremely strong 
(gi — > 00), the JT distortion is large, higher order terms (neglected in (P) will in prin- 
ciple become relevant. In such a case, the final state tends to approach closely a static 
JT symmetry-broken distortion. When the coupling is instead large but finite, there is a 
semiclassical region, where quantum fluctuations are able to connect classical equivalent 
valleys via tunneling, and so restore the full symmetry while also generating characteristic 
low-energy tunneling excitations in the spectrum. The intermediate coupling region gi ~ 1 
is, as usual, the less predictable. However, the weak coupling limit (gi — > 0) finally simplifies 
again into a perturbative problem. If either g 4 or g 2 alternatively vanish, we will have the 
special cases of © T)^ or T>^ © T)^ JT respectively. A special case we shall consider 
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in greatest detail below is that of equal coupling of degenerate oscillators u 2 = ^4, gi = #4- 
This paper is organized as follows. The weak-coupling region g 2 , g$ « 1 is discussed 
in Sect. [H], obtaining approximate analytic expressions for the splittings and shifts in the 
spectrum. Next, we shall move to the opposite, strong-coupling case. Here, the classical 
JT minima of the Born-Oppenheimer potential energy surfaces are studied first (Sect. pTT] ); 



then, in Sect. [TV], the ionic kinetic energy is introduced semiclassically to obtain a descrip- 
tion of the low-energy quantum excitations in the strong-coupling limit. Exact numerical 
diagonalizations will be used to shed light on a few particular features of the intermediate- 
coupling regime. Finally, the relevance of these preliminary theoretical results to the physics 
of fullerene cations is discussed in Sect. |V|. 



II. THE WEAK-COUPLING LIMIT: PERTURB ATIVE APPROACH 

The weak-coupling limit in a JT problem is a coupled situation in which both the elec- 
tronic and the vibronic degrees of freedom need to be treated on equal footing, as equally 
"fast". Their cooperative dynamics is strictly regulated by quantum mechanics, and no 
"classical degree of freedom" can be singled out in this limit, as opposed to the Berry-phase 
semiclassical approach. Nonetheless, the limit of small g 2 , g± is easily treated by perturbation 
theory, and the results are useful in providing an exact match to the semiclassical spectrum 
on the opposite side of the range of coupling. 

The JT coupling (Q) splits and shifts the harmonic levels to second-order in g^. We 
compute here these shifts for both the fivefold-degenerate ground state and the 70-fold 
degenerate one-vibron excitation in the 1-hole case and for equal frequencies Tiuj 2 = Uoj& = 1. 
As we have two small parameters, g 2 and g^, and two different perturbing Hamiltonians, 
Hff_ v and H^} v , there is an infinite number of choices of the small perturbing parameter, 
depending for example on the ratio g 2 /gA- We compute the energy shifts for three special 
cases: g^ = 0, i.e. perturbing with respect to H^} v , g 2 = 0, using H^_ v alone, and finally 
#2 = #4 = g, the equal coupling case. The coefficients for the three cases are collected in 
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Tab. |I|, along with the appropriate SO (3) symmetry labels. The energy of a state (in units 
of Too) is thus given by 

(number of vibrons) + (coefficient in table) x g 2 + 0(g 4 ) , (5) 

where the number of vibrons is one for all the states except for the ground state at the first 
line. 

We note on the second column five states, in the 1-vibron multiplet, whose shifts are 
— \g\, he. the same as for the ground state. These states are readily identified as the 45- fold 
degenerate state corresponding to the T>^ vibron fundamental state (here uncoupled, since 
#4 = 0). Parallel considerations apply to the uncoupled T>^ vibron in the third column, 
recognizable by the shifts of — Jj^J. 

The equal-coupling case is interesting from the point of view of symmetry. Pooler studied 
the Hamiltonian (|j) from an algebraic point of view |S| , discovering that the equal frequen- 



cies - equal coupling case uj 2 = uj^, g 2 = g±, is actually characterized not only by the obvious 
SO (3) symmetry, but also by rotation symmetry in a 5-dimensional space (SO (5) group). 
In the language of the SO (5) irreducible representations, characterized by two quantum 
numbers [l,m] ||24|] , the T>^ hole is classified as a [1,0] irreducible representation of SO (5), 
while the 14- fold degenerate vibration [V^ ®T>^) is a realization of the [2,0] representation 
of SO (5) [p5[| . In the fourth column of Table ||, we note that the one- vibron level in the 
equal-coupling case splits into three levels with degeneracies 30, 35 and 5 respectively. The 
levels of degeneracy 5 and 30 are readily seen as representations [1,0] and [3,0] of SO (5) 



|. The 35-fold degenerate level, by table 10-3 of reference p4| is identified with a [1,2] 
representation, as indicated in the last column of Table |I[ As this large symmetry is an 
exact feature of the Hamiltonian, it is not limited to the perturbative regime, but will be 
carried along to all values of (equal-) coupling, and to the semiclassical limit (Sect. [TV]) in 
particular. 

To conclude, we note that, while the L = 0, 1, 5, 6 levels are exactly the same in the 
three perturbative regimes summarized in Tab. J|, the L = 2,3,4 levels (marked in Tab. 
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II] with symbols v,*), in the three cases, are instead different linear combinations of the 
same states: this fact has a very simple explanation. The L — 0, 1 levels uniquely come 
from the coupling of the T>^ hole to the T>^ vibron, and the L = 5, 6 levels derive from 
the coupling of the T>^ hole to the T>^ vibron. The L = 2, 3, 4 levels can come from either 
of the previous couplings, and there is no a priori reason why they should not mix among 
themselves. 



III. TOPOLOGY OF THE JAHN-TELLER MANIFOLD AND THE BERRY 

PHASES 

In this section we want to understand if the JT hole-vibron coupling in Ceo positive ions 
possesses a Berry phase or not. In order to do that, the canonical steps are the following: (i) 
find the manifold, in the space of classical ionic coordinates, formed by the set of equivalent 
minima of Hamiltonian (^) (the minima of the Born-Oppenheimer (BO) energy surface); (ii) 
study the adiabatic transport of the electronic ground state wave-function - in this case a 
unit vector spanning a five-dimensional sphere - while the ionic coordinates are taken along 
a closed path on their manifold. Representation (gj) of Hamiltonian (H) is inconvenient for 
this purpose, and we first perform a transformation of the electronic and vibronic operators 
to a suitable representation. 

A. The real representation 

We start defining a vector of hole-destruction operators, 

a = ( (6) 

In terms of this vector, the interaction Hamiltonian can be written as 

H&^CI-Bl.C, (7) 
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where indicates the adjoint transpose of (|6]), and B 2 and B4 are 5x5 matrices whose 
elements are combinations of vibrational coordinates. Next, we apply the following trans- 
formation to the electronic operators 



Clcr 



C„ 



' % i ^ 



1 

72 



1 -1 

y/2i 

-i i 

1 0-10 



(8) 



C2fr 
C3<t 
CAa 

K ~c 5aJ 

The transformed vector C a contains the phase factors necessary to give a suitable symmetry 
to the transformed 5-matrices. 

The final step converts the vibron operators into their real coordinate representation: 

L 



q L ,m= E M m , A (^ iA +(-l) A 6 Ll - 



X=-L 



where 



(9) 



M m ^o = (2sign(m)) 2 (<5 m>A + (-l) m sign(m)5 m _ A ) 

and the masses of both oscillators are taken as unity 

We rewrite H and Hh_ v after the final transformation as 

Ho = E =r E H2,. + 



L=2,4 " l=-L 



where the matrices B 2 and B4, the transformed of B 2 and B± under 
of qL, m (|) read 



(10) 



(11) 
(12) 

expressed in terms 
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B 2 = -= 



1 

7i 



292,0 





292,2 


\/3g 2 _i 


-V3?2,l 





2(79 n 


2(79 9 

Z, — Z 


— VO(79 1 


— VO(79 1 


292,2 


2 92,-2 


-2g 2 ,o 


"92,-1 


-<?2,1 


\/3g 2 -i 


-V3g 2 ,i 


"92,-1 


-g 2 ,o + v / 3g 2 , 2 




-^92,1 


-\/3<? 2 -i 


-92,1 


-v / 3g 2 ,- 2 


-92,0 - V3q 



(13) 



and 
73 4 = 



56 



5 94,o + Vl4g 4 ,4 

^94,-4 
^94,2 

g 4 ,-i + V7q±-3 



5 04,0 



14g 4 _ 4 ^94,2 
V / 14g4,4 V / 6g4,-2 



i/4,-2 



^4,-1 + V7q4,~3 
-94,-1 - ^94,3 
5 94,0 \/l2<74 _i 



"94,1 



\/7g4,3 a/1294,-i -^94,0-^94,2 



y -94,1 + ^94,3 -94,-1 + ^94,-3 V / 1294,l 



4,-2 



-94,-1 + ^94,3 
-94,-1 + ^94,-3 
^94,1 

VSq^-2 
■yjf 94,0 + \/894,: 



(14) 

\ 



B. The Born-Oppenheimer energy minima 

The classical (purely static) solution of the JT problem consists in finding the minimum 
of the Born-Oppenheimer (BO) potential energy defined as the sum of the potential term 
in (|TID plus the lowest eigenvalue of the hole-vibron interaction matrix 



M (92, 94) = 92— B 2 + g±— B A 



(15) 



in the (5+9)-dimensional space of the classical qL, m coordinates. This classical minimum is 
generally not unique. More typically, there will be continuous sets of minima, forming well 
defined manifolds embedded in the multi-dimensional coordinates space ||27|| . We call one 
such set of minima of the adiabatic potential Jahn- Teller Manifold (JTM). The topology of 
this manifold crucially affects the possible presence of a Berry phase. We study the JTM 
following the approach of Opik and Pryce more recently reformulated by Ceulemans in 
a clear and useful form ^Z7\ . 
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Consider a generic one-fermion state 

5 



I*7> = E«,tI°>> ( 16 ) 

where without loss of generality, the 77^- can be taken real, with the normalization constraint 

X>? = 1- (17) 
i=i 

The quantum amplitudes ff live therefore on a S4 sphere in 5-dimensions, which we indicate 
as the electronic sphere (ES). 

Within this notation, it is convenient to rewrite the classical potential energy as 

hi > L 

v(g 2 , v) = E D q h + *f • &) • »f (is) 

L=2,4 ^ i=-L 

at fixed distortion <f 2 , and electron state ff. Minimization of V with respect to ff yields, by 
definition, the (adiabatic) BO potential. If we keep ff fixed instead, and minimize flT8|) with 
respect to q — (q 2 , q\), then we get the extremal points of the BO potential in configuration 
space as 

q=q(ffi (19) 

Inserting ([19|) in (0), we obtain a function V of the electronic coordinates alone, hence 
a function defined on the ES, not on the space of distortion like the BO potential. We still 
have to ensure that the electronic vector is an eigenvector of the interaction Hamiltonian. 



Ceulemans [|27j] realized that this condition is equivalent to a second minimization procedure, 
on the ES; moreover, he recognized that the V function is isoextremal to the BO surface, 
that is, its extremal points on the ES correspond (in a 2:1 way) to extremal points of the 
BO surface. Then the JTM has the same symmetry properties of the set of minima of V. 

We apply this procedure of reverse-order minimization to the present system. We list 
here the distortions q 2 and 04 which minimize V at fixed if on the ES: 

2 1 [E 

q2,2 = -^92vm - 2 y 7^2(^4 - vl) 
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3 / 1 

92,1 = y ^92{V2V4 + ViV5) + ^92VsV5 

g 2 ,o = —JfM + V2 - vl) + ^92(vl + vl) (20) 

92,-1 = \j^g2(V2V5 ~ Vm) + -^92V3V4 

2 [3 

92,-2 = ~^92V2V3 + y y92» 

94,4 = —giivl - vi) 

94,3 = -^9i{WA - V1V5) 

94,2 = -^gmm + -^vMl - vt) 
94,1 = y yi9AV2Va + vivs) - y jV3Vs 

' 9a{Vi-4- Qvl + H + H) (21) 



1.0 



2V35 



94,-1 = -y ^9i(vm + V2V5) - J^vm 

[3 2 
94,-2 = -\j^94V2Vs + ^gmvs 

94,-3 = -^g^vm + V2V5) 
94,-4 = -g4ViV2 

We can choose a (hyper-polar) representation for the electronic coordinates such that the 
unitarity condition ( |17|) is automatically satisfied 

r/i = sin 6\ sin 2 sin 3 cos 04, 
r] 2 = sin 6 1 ! sin 2 sin 3 sin 4 , 

7/3 = cos 0i, (22) 

774 = sin0! cos0 2 , 

775 = sin 0! sin 2 cos 3 . 

Here, 0; i=l,2,3 are the azimuthal angles and 04 the anomaly. Substituting (^) in ( |2"0D and 



(pi|), and then back into ( J1S[) we obtain an expression for the minimum .Ej^ of V: 
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E JT = - — (5^2 + 9gl^4)- 



(23) 



and for the JTM coordinates 



1 3 , „, 

m\ = -^92 , m = -^9i ■ (24) 



C. Is there a Berry phase? 

As Eq. ( p3|) shows, Ejt is independent of the r\ parameters, and is therefore a constant on 
the whole ES. This is a nontrivial consequence of the large amount of geometrical freedom 
which is allowed to the molecule by the large degeneracy of the vibron modes: on the 
contrary, when the dimension of g-space is less than that of the ES (like for example in 
the h <S> G JT system - relevant for the coupling of h u with the G g (3) mode, which has 
no H g partner), Ejt is not flat, but is characterized by some amount of corrugation. The 
isoextremality property of the T>^ ® (T)^ © T>( 4 )) implies that the JTM, must have locally 
the same structure as the ES, here S 4 , a sphere in five dimensions. 

However, upon closer scrutiny, we note a topological difference between the ES and the 
JTM, due to the 2:1 correspondence mentioned above. In particular, the ES is a projective 
space, as opposite points on the sphere correspond to the same point in the g-coordinates 
space (they belong to the same "ray"). This is equivalent to the fact that normalized 
eigenvectors of a symmetric matrix are defined up to a sign. The JTM does not share this 
property, since there are no pairs of points on the JTM that generate the same fj. Here is 
the source of the Berry phase: to each path in the JTM there corresponds a path on the ES, 
but the two are topologically different. In particular, a loop in the JTM may correspond 
either to a closed path, or to a path from a point to its opposite (antipode) on the ES. In 
general, the latter path cannot be reduced to a point, and there will be a topological effect 
(Berry phase). This fact is illustrated in Fig. [1] on an ordinary S2 sphere, in 3-dimensional 
space, representing the ES. The line V corresponds to a closed loop on the JTM: there is no 
way to make point B coincide with point A in such a way that then we can just shrink the 
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path and have it coincide to a trivial one such as T±. However, we are going to encounter 
situations where apparently nontrivial loops become shrinkable to a point via singularities 
in the mapping between the ES and the JTM. 

Of the most general set of coupling parameters, we consider three special cases: (a) 
® £>(2) ; i.e. g 2 0, 04 = 0; (b) <g> i.e. g 2 = 0, g 4 ^ 0; and finally (c) 

P (2) <g> (X> (2) © P (4) ), i.e. cj 2 = u A , g 2 =g A ^ 0. 

(a) "D^ ® T>( 2 \ Here, the <f — ff mapping is given by Eq. (|20"1) and (p2|) . A class of closed 
paths on g-manifold does, as always, transport ff into its antipode — ff. However, we find 
in this case that a whole ("equatorial") line joining two antipodes on the ES, defined by 
9\ = 9 2 = 9 3 = vr/2 (i.e. r/f + r/| = 1, r/ 3 = 774 = 775 = 0) and parameterized by 6*4 = — > n, 
in fact corresponds to a single point q 2 = (0,0,-^,0,0) in g-space. The antipodes are 
therefore topologically joined, and there are no loops on the JTM that cannot be deformed 
to a point. Hence, there is no Berry phase in this case. 

(b) T>^ © T>( 4 ). The q — ff mapping is now through Eq. (|2T| ) and (p2|). A closed loop 
on the JTM transports ff into —ff if it encircles a conical intersection in configuration space. 
This is clearly the case of the line described in (a), here corresponding to 

<?4,4 = -7^4C0S(26) 4 ) 

94,o = -— t=#4Cos(20 4 ) (25) 
g 4 _ 4 = -^ 4 sin(20 4 ) 

Now only the q antipodes are topologically equivalent, while the loops that join them are 
nontrivial: there is a Berry phase of ir in the electronic transport in this case. 

(c) T>( 2 ) © (X>( 2 ) © X>( 4 )). In this case the q — ff mapping passes through both Eq. ( |20|) 
and Eq. (PH). The Berry phase of case (b) is therefore conserved upon addition of the T>^ 
mode. This case has an instructive analogy to simpler JT-coupled systems. 
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D. "Propagation" of Berry phases 



In case (c) above, we add to a Berry-phase entangled DJT system an additional L = 2 
vibron mode, also coupled to the same electronic state. It is interesting to consider what 
happens to the Berry phase when the added mode, on its own, carries no Berry phase. The 
enlargement of the <f-space introduces a large freedom for creating new loops in q. One might 
in principle think that some of these loops could manage to avoid the conical intersections 
of BO levels, thus healing the connectivity properties of the JTM, and short-circuiting out 
the Berry phase. We found instead that this is not the case, and T>^ ® (T)^ ®D^) retains 
the Berry phase of £> (2) ®V {i K 

An example of analogous situation is the t®(T@E) problem, extensively studied in the 
literature []n|[n|: the t®T problem is Berry-phase-entangled, while the t <E> E is not. It is 
well known that t®(T@E) problem retains the same topological properties, including the 
Berry phase of 7r, of the t <g)T system, the effect of the E mode consisting essentially only 
in flattening the JTM, and removing inter-minima barriers. 

The reason why new loops seem to be ineffective in destroying the Berry phase in the 
total vibron space is that this space is the cartesian product of the single-mode spaces: the 
degeneracy points, opening "holes" in the space of the first vibron, propagate (locally, at 
least) in the global space along the second-vibron directions. The global conservation of 
connectivity in the whole product space is a nontrivial topological problem, whose general 
outcome is not known to us. In both examples above, we simply note that things go as if 
the Berry phase carried by the first mode is strictly conserved. 

In the present situation, the coupled T>^ vibrational mode does not present any topo- 
logical effect (a), whereas the as seen in case (b), carries a Berry phase. Thus the 



overall T>^ cg> [T>^ © T)^J system is affected by topological effects, and this remains true 
even for different couplings gi 7^ g& 7^ 0. 

The fullerene cations provide a third example of this behavior. The eg) T>^ is nat- 
urally resolved in its icosahedral composition h® G plus h ® H. In the light of the above 
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considerations, we expect that in this system the Berry phase should be carried by the G 
mode ||28|| . That this is indeed the case, can thus be seen as a natural consequence of the 
absence of Berry phases in h H (i.e. T>^ T>^), and its presence in 

IV. SEMICLASSICAL QUANTIZATION AND THE EFFECTS OF THE BERRY 

PHASE 

So far we have treated the ions as classical. We study here the consequences of the 
discussed symmetries of the JTM and of the presence/absence of the Berry phase on the 
quantum spectrum of the coupled system, by quantizing the semiclassical ionic motion. The 
general form of the ion kinetic energy is 

L=2,4 Z l=-L U( 1L,1 

It is convenient to use the classical language for the derivation of the kinetic operator re- 
stricted to the semiclassical region: 

1 L=2,4 U L l=-L 

In the semiclassical limit, valid when gi are very large, the ionic motion is confined to the 
JTM, itself given by Eq. fl24l) . Hence, only the tangential part of the kinetic energy (|26D is 
relevant, which, upon introduction of Q = ^ — 1 can be written 

K = ^-E^ G % > ( 28 ) 
where G'- 7 is the metric induced on the JTM by the electronic parametrization: 

G» = £ (1 + 5 LA Q £ deSLAi^i ■ (29) 

L=2,4 l=-L 

Not much progress is possible in this form, due to the general low-symmetry of the kinetic 
term. However, in the special case of equal frequencies {0J2 = oj^ = u: ( — 0) and equal 
couplings, where the true symmetry of the problem is SO(5) p3|, the metric is rewritten as: 
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G*' = d di qde 3 q (30) 
Thus, the kinetic energy corresponds to that of a free particle on a sphere in 5 dimensions: 
K = I—{6i + sinO'/ + sin0i 2 sin0 2 2 3 2 + sin 0i 2 sin 2 2 sin 9 3 % 2 ) . (31) 
I is the 'inertial momentum' of the system 

/ = l?1 2 = \a 2 (32) 

The corresponding Laplace-Beltrami operator (describing the free motion of a quantum 



particle on the manifold parametrized by the 9% variables) is given by the formula [29 



A = -l=d et G^VGd ej , (33) 
V G 

where G is the absolute value of the determinant of the metric tensor. This operator is 
the usual representation of the quadratic Casimir operator of SO (5) on the coset manifold 
SO (5) /SO (4), which is actually the sphere S 4 : this is a consequence of the exact SO (5) 
symmetry of the equal-coupling case. 

The particle moves freely on the JTM, with harmonic oscillations orthogonal to it. The 
corresponding energy levels are therefore 

h 2 10 1 

E = E JT + —1(1 + 3) + tuv J2 K + ^) (34) 

^ 7=1 ^ 

Ejt is the static JT energy (^). The last term is due to the vibrations orthogonal to the 
JTM. The 1(1 + 3) term corresponds to the eigenvalues of the Laplace-Beltrami operator 
on the sphere; I is the label of the angular hyperharmonics on S* 4 , that realize the [1,0] 
representations of SO (5), whose degeneracy is given by 

= (2l + m + m + ll 

6 

The Berry phase imposes precise selection rules on these levels: since the one-hole state 
changes sign after a closed loop around the JTM, and the overall wave function needs to 
be single- valued, also the vibrational wave function has to change sign. Odd-/ states have 



this properties, and these are the allowed states that survive the Berry phase selection rule. 
This implies that the ground state is I — 1 (fivefold degenerate) and the first excited state 



is I = 3 (30-fold degenerate). Results of exact diagonalization on a truncated basis P,|30|JTl 



* = E a oA |0> + a ia b\cl |0> + a ija b\b]ci |0> + ... (36) 

shown in Fig. ^], confirm the spectral structure described here for the semiclassic limit. At 
all values of the coupling g, as expected, the degeneracies are those characteristic of SO(5). 
We have done also numerical diagonalizations in the more general case of different coupling 
and/or frequencies: the spectra, as expected, display lower - SO(3) - symmetry, but a 
qualitatively similar behavior. 



A. Two holes in spin-singlet state 

In the spin-singlet n=2 holes case, 1 Cg f we find two main differences with the n=l 2 Cq" 
case. First, the average distortion \q\ = 2\f^g is twice as large as in the n=l case, / and 
the Ejt energy being 4 times as large. Secondly, and most importantly, there is a change in 
the selection rules. The two holes live in the same single- hole eigenvector of 7i(<f2 ; $i)> the 
electronic wave function being the product of two equivalent eigenfunctions for the space 
coordinates and an antisymmetric part for the spin degrees of freedom. Since each hole 
state changes sign after a closed loop in the JTM, the overall hole state does not. Thus the 
vibrational wave function must also be single- valued in order to have an overall single- valued 
wave function for the system. As a consequence, only even-Z states survive: in particular, 
the ground state is / = (non-degenerate) and the first excited state is / = 2 (14-fold 
degenerate). Exact diagonalization results, shown in Fig. |2](b), confirm that the symmetry 
and, to a certain degree the spectral structure described here for the semiclassic limit are 
retained in the weak-coupling limit. 
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B. The 2 (g> (2 © 4) unequal couplings case 



In the more general non-equal coupling case g 2 ^ #4, including also the g 2 = case, 
we could try to use the metric as for the equal coupling in order to obtain an analogous 
expression for the kinetic energy. Unfortunately, when the couplings are different, we do not 
get a treatable expression like (|3T|) , and we cannot follow the same route for quantizing the 
system. 

We recall instead that the structure of the JTM is the same as that of ES, and in 
particular that SO (5) symmetry is retained. The expression of the minimum potential 
energy (^3|) suggests a flat JT trough, for any values of the frequencies and couplings. 
Therefore, it is likely that the free motion of the particle on the JTM is described again by 
the quadratic Casimir operator of SO (5), retaining for the low-energy excitations a spectral 



structure similar to Eq. (p4|) . 

However, although the full SO (5) symmetry is retained in the g — > oo semiclassical 
limit, for any non-infinite value of the couplings, the SO(5) representations will be split 
into representations of the SO (3) group R5I| |. From this point of view we can see the dif- 
ferences between the frequencies and between the coupling strengths as symmetry-reducing 
perturbations, that become effective when moving away from the strong-coupling limit. 



C. The 2 2 case 

This case, corresponding to h u ® H g in a Cgo cation, is very interesting because of the 
peculiar features described above, i.e. absence of topological singularities and, in particular, 
absence of a Berry phase. Furthermore, this case is worth additional analysis, since it 
provides an alternate limit of our approach for the different coupling case (of which this is 
an extreme case). 

We treat again the system as a quantization problem on a five-dimensional sphere. In 
this case, it is possible to find the extremal eigenvalues of the interaction matrix 7i(q2) = 
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|g 2 ^2-£>2, as solutions of the characteristic polynomial equation for B 2 : 

A 5 - |g 2 | 2 A 3 + 7/ (g 2 )A 2 + ^|g 2 | 4 A - 4/ (g 2 )|g 2 | 2 = (37) 



49' 



where 



/<)(&) = ~^/^^ q2 ^ q l 1 + g 2,~l) ~ 6 ?2,0(g 2 ,2 + ?2 -2) 

+2g| + 3v / 3g 2 , 2 (g 2 ,i - 92,-i) ~ 6^2,1^2 ,-192,-2] • (38) 



This polynomial can be decomposed as 



(A 2 -^ 2 | 2 ) 



The extremal eigenvalues of B 2 are 



A s -?|fe| 2 A + 7/o(fe) 



(39) 



2 

A1/5 = T^|g 2 | , (40) 

in accord with the SO (5) symmetry of the JTM in the 5-dimensional T>^ space. The three 
remaining intermediate eigenvalues can be expressed as 

2 2tt\ 
A 2 = N 7f cos( v -- Y j 

A 3 =|g 2 |- 7f cos( v - 

A 4 = |g 2 |^=cos ( - + — ) (41) 

where 



y/7 V3 3 y 



cos 



* = -49V7^ . (42) 



2|g 2 | 

In order to obtain the eigenenergies of Ti,(q 2 ) one has to multiply these eigenvalues by a 
factor \g 2 liuj 2 . By inspection, A 2i 3 5 4 are bound in the [Ai, A5] interval. However, it is possible 
to find many points on the JTM such that Ai = A 2 , for example the "south pole" qo/\q\ = —1 
and a one-dimensional manifold on the "equator" q = 0. Since Ai is constant and A 2 is 
analytical in q 2 , in particular around the tangency points where the two eigenvalues get 
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degenerate, these degeneracies are not conical intersections of the two adiabatic surfaces. 
The two surfaces have therefore just a contact of second order at these points. The absence 
of conical intersections implies the absence of any topological effect 0. 

Because of the absence of a Berry phase, already discussed in Sect. |TTTj, the strong coupling 
behavior of the system, described in terms of the SO (5) representations, is such that we 
expect a [0,0] (non degenerate) ground state, a [1,0] (fivefold degenerate) first excited state, 
a [2,0] (14-fold degenerate) second excited state, and so on. Moreover, absence of topological 
effects makes the two-holes case not essentially different from the one-hole problem in the 
strong coupling limit. Of course, the static JT energy and inertial momentum is just four 
times larger than in the one hole case, but the level sequence in the strong-coupling limit is 
the same. 

Coming back to the single-hole case, let us now move away from the strong-coupling 
limit, to discuss the more general case of an arbitrary coupling g 2 . We might naively expect, 
at first, the ground state symmetry to be preserved for arbitrary g 2 , as was the case in other 
DJT problems examined previously ||. However, a closer inspection reveals that this is not 
true, and rather that there must be a level crossing as a function of g 2 . A first indication 
comes from the fact that the ground state must in fact be L = 2, i.e. fivefold degenerate, as 
g 2 —* 0, reflecting just the uncoupled hole (h u ) degeneracy. A second hint comes from the 
tangency of the lowest with the second BO surface, noted above. Could this tangency have 
a role in the expected level crossing from L = 2 at g — > to L = at g ^ oo? To study this 
possibility, we resort to numerical diagonalization of Hamiltonian (|3]) on a truncated basis 
(|36|), for finite values of g 2 . As pointed out above, for finite coupling the symmetry of the 
system is strictly SO (3). Therefore, we expect the large degeneracies characteristic of SO (5) 
symmetry to build up from converging SO(3) states for large g 2 . The [0,0] and [1,0] states 
map directly on the T>^ and SO (3) representations respectively ||31|| , while the SO (3) 



composition of [2,0] is © V^. Therefore, we seek the convergence of two such levels in 
the strong-coupling limit. 

The numerical solution of the 1-hole problem (see Fig. |||(a)), shows indeed the correct 
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order of the levels for large couplings (g 2 ^ 8), whereas for g 2 < ~8 the [0,0] and [1,0] 
levels cross and become inverted. This level inversion can in fact be precisely seen as a 
consequence of the tangency of the two adiabatic potential surfaces described above. The 
two BO sheets lie close in a significant region of configuration space. Treating the upper 
sheet as a perturbation for the eigenstates of the lower sheet, the energy of all these states 
is pushed up by this correction. The / = (isotropic) states feel more efficiently this 
perturbation than the / > hyperspherical harmonics, since they can arrange their nodes 
so as to avoid the tangency line. In conclusion, the effect of the tangency is to increase the 
energy of the I = level more than that of / = 1, so that, for coupling not too large, g ^ 8, 
the latter is the ground state. 

In the 2-holes (spin singlet) case (Fig. 0(b)) this inversion is not present, the ground state 
being nondegenerate (T>(°>) for any coupling g 2 . Again, from a weak-coupling approach, this 
is a simple consequence of the order of splitting of the d 2 -fermion spin-singlet states upon 
perturbative JT coupling with a T>^ vibration, that is (from lowest to highest) D", T>^\ 
p( 2 ). From the strong-coupling side, the only observed crossover happens between the T>^ 
and T>^ low excitations above the ground-state. This means that the BO picture of the 
two holes in the same single-hole orbital, although appropriate to the strong-coupling limit, 
is rather poor for intermediate JT coupling, when the two fermions manage to arrange in 
a cooperative state preserving the global T>^ symmetry despite the tangency of the lowest 
BO sheets. 



D. The h u <8> G g case 

Looking at Table | we can see that there is a single G g left, that does not pair with 
any H g mode. We must therefore treat the coupling of mode G g (3) according to the true 
icosahedral picture. Actually the static h®G JT problem has been treated in the literature 
. The set of absolute minima of the system if made up of ten isolated points in the four- 



dimensional G g space. Tunnel splitting among these ten BO valleys give rise to vibronic 
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states of symmetry A, G, and H. In Ref. [28fl the diagonalization of the corresponding 



matrix has been performed including the proper provision for the Berry phase. The result is 
that the H quintet is the true ground state, followed by the G quartet and the A singlet in 
the order. Thus, the Berry phase again acts in such a way to enforce a fivefold degenerate 
ground state all the way from weak to strong coupling, in analogy with the h © (G © H) JT 



problem, and in contrast to the h® H case discussed in Sect. [IV C| . This confirms that, as 



anticipated, the Berry phase in the h <S> {G © H) JT problem comes from h®G sector. 

V. DISCUSSION 

We have studied the symmetry aspects and the possible presence of a Berry phase in Cg 
and Cgo" in the linear and spherical approximation. We find that while 1 Cqq has no Berry 
phase, 2 Cg does, and is very rich and intriguing. 

Clearly, the present analysis is only a first step. In reality, fullerene is not a sphere, but 
a discrete icosahedron, and the modes are split into H g © G g modes. Nonetheless the 
real shape of the molecule can be seen as a perturbation on the spherical picture (Tab. |): 
the SO(5) and SO(3) levels are split but the selection rules are still present. The Opik- 
Pryce approach can still be applied in order to obtain the JTM in icosahedral symmetry 
|32]| . Quite generally the BO potential will warp in such a way to originate a discrete set of 



minima, instead of a continuous 4-dimensional manifold. The molecular system will execute 



tunneling among these minima. However, as shown by Ham [33[] , the Berry phase imposes the 
degeneracy of the ground state (irrespective of the presence and height of barriers between 
the BO minima), which still is fivefold degenerate {T>^ — > H). Therefore, we expect the 
ground state symmetry to be 2 H U for real CgQ, which has a Berry phase, and l A g for singlet 
Cgo , which does not. 

A quantitative study of the properties of Cg is clearly beyond the scope of the present 
work. In order to accomplish that, a full set of e-v couplings gu g {i) i = 1, -..8 and gc g (i) i — 
1, ...6 will have to be known (either from calculations, or fit to experiments), and inserted 
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into a diagonalization procedure similar to that previously carried out in the case of Cqq" 
|p| j30|JTT|1 . While that remains to be done, we can list, for the time being, some tentative 
physical implications of our results. 

(i) Low-energy electron detachment. In a way analogous to that proven for the case of 
Cg , one can surmize that Cg~ cannot have any A g low-energy states. It follows that if an 
electron is detached from Ceo into a low-energy state, the s-wave cross-section will be zero, 
and there will be a centrifugal barrier in detachment too, similar in nature to that pointed 
out for attachment [|l(|. This effect could be related to the thermoionic anomalies pointed 
out by Yeretzian et al. |34| 

(ii) Low- energy tunneling excitation. In the spherical approximation, the lowest excita- 
tion of Ceo is a 30-fold degenerate SO(5) "/ = 3" tunneling multiplet, as required by the 
Berry phase. Split, as it will be by the lower SO(3) symmetry, as well as by the true 
effects, this multiplicity might still be identifiable spectroscopically. 

(iii) photoemission spectrum. In the valence photoemission spectrum of Cgo — * 
Cqq + the electron kinetic energy distribution contains very direct information - in the 
so-called final state approximation - about vibron and hole-vibron coupling in Cg~ . The 
outgoing electron shakes up multi- vibron excitations of C^ , with a relative probability which 
can be calculated in principle, in a manner similar to that used by Gunnarsson |3(| for Cg . 
Very recent data [35] indicate that these shake-ups do exist and are in fact stronger for Ceo 
than for Cg . 

We hope to return to these topics in subsequent work. 
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FIGURES 

FIG. 1. A three-dimensional representation of the (truly five-dimensional) electronic sphere 
(ES): the path T connects opposite points (antipodes), which correspond to the same point on the 
JTM. Thus r maps, through Eqs. ( p0| ) and (|2l|), on a closed loop in the JTM. This class of loops 
is usually entangled with a nonzero Berry phase, while trivial loops, such as T± can be smoothly 
deformed to one point, with associated zero Berry phase. However, in the case 54 = 0, due to the 
singularity of mapping (pop, a whole equatorial line A maps to a single g-point. In this special 
case, point B identifies with point A on the JTM, thus T becomes equivalent to T±, and no Berry 
phase is possible. 

FIG. 2. The low-lying states of the spectrum of the L = 2 ® (L = 2 L = 4) ([l,0]<g>[2,0], 
in SO (5) notation) JT system, as a function of the coupling g = g\ = g\- Symmetry labels are 
indicated according to the SO(5) group representations. The truncated basis includes up to 7 
vibrons. Panel (a): the levels originated from the 0- and 1-vibron states, for n = 1 holes (see also 
Table || for the small-g limit). Panel (b): the three levels originating from the 0- vibrons multiplet 
for n = 2 holes. The tenfold-degenerate [0,2] orbital state is in a triplet spin configuration. 

FIG. 3. The low- lying states of the spectrum of the L = 2 L = 2 coupled JT system as a 
function of the coupling g\. We include in the diagonalization states with up to 35 vibrons. For 
n = 1 fermions included (panel a), the reference level is taken as the lowest L = 2 state, that is also 
the ground state for 52 <~ 8. For n = 2 fermions included (panel b) the energies are excitations 
above the L = ground state. Only spin-singlet states are plotted. 
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TABLES 


S0(3) 
L 




Calculated 
cm -1 


22 


] Calculated 36 
cm -1 




2 


H 9 (l) 
H 9 {±) 


214 
727 




261 
775 




3 


G g (3) 


557 




594 






H 9 (2) 


387 




435 




4 


Gg(l) 


374 




482 




Hg(6) 

G B (4) 


1180 
774 




1208 
1047 




5 




1091 




1098 






H 9 (3) 


516 




730 




p. 
u 


G g (3) 


623 




781 




H 9 (7) 
G g (5) 


1477 
1412 




1394 
1314 




7 


H 9 (S) 
G g (6) 


1617 
1694 




1573 
1479 





TABLE I. Experimental energies of the JT-active G g and H g modes, organized according to 
their spherical parentage p^| . The numbers in parenthesis in the second column are the position 
of the mode in order of increasing energies. 
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S0(3) 
L 

2 (GS) 


3t 
4* 
6 
1 

2* 
3t 
4* 
5 
2* 



AE 



ii 

'28 



_9_ 

'28 



+ 



56 



AE 
~9i 



__9_ 

20 

__9_ 

20 

121 
280 

141 

'280 

_X 
10 

__9_ 

20 

__9_ 

20 

__9_ 

20 

__9_ 

20 

13 
40 

1 

280 



AE 



7_ 

' 10 



> — 



19 
20 



SO(5) 

Ml 

[1,0] 



[3,0] 



23 
40 



[1,2] 



_3_ 

'40 



[1,0] 



TABLE II. The coefficients of shifts of the vibronic levels for second-order perturbation in the 
couplings (72 and 54, for the 0- and 1-vibron multiplets of levels, for n = 1 holes. In the first column 
we list the labels of the states according to spherical symmetry; the second and third column report 
the results of applying as a perturbation H^} v and H^} v alone respectively; in the next column 
the perturbation is done with equal couplings Q2 = 9a = 9- In this last case the true symmetry of 
the system is SO(5), that we evidence with the labels in the last column. The states marked with 
t and are multiply present in the 1-vibron multiplet, and therefore admix differently with the 
three different perturbations. Results in units of Hu)2 = hoj^. 
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